Draft version October 31, 2011 

Preprint typeset using I^T^}X style cmulatcapj v. 11/10/09 



ON THE LACK OF TIME DILATION SIGNATURES IN GAMMA-RAY BURST LIGHT CURVES 

Daniel Kocevski ^ and Vahe Petrosian ^ 

Draft version October 31, 2011 

ABSTRACT 

We examine the effects of time dilation on the temporal profiles of gamma-ray burst (GRB) pulses. 
By using prescriptions for the shape and evolution of prompt gamma-ray spectra, we can generate 
a simulated population of single pulsed GRBs at a variety of redshifts and observe how their light 
curves would appear to a gamma-ray detector here on Earth. We find that the observer frame 
duration of individual pulses does not increase as a function of redshift as one would expect from 
the cosmological expansion of a Friedman-Lemaitre-Robertson- Walker Universe. In fact, the duration 
of individual pulses is seen to decrease as their signal-to- noise decreases with increasing redshift, as 
only the brightest portion of a high redshift GRB's light curve is accessible to the detector. The 
results of our simulation are consistent with the fact that a systematic broadening of GRB durations 
as a function of redshift has not materialized in either the Swift or Fermi detected GRBs with 
known redshift. We show that this fundamental duration bias implies that the measured durations 
and associated E^so estimates for GRBs detected near an instrument's detection threshold should be 
considered lower limits to their true values. We conclude by predicting that the average peak-to-peak 
time for a large number of multi-pulsed GRBs as a function of redshift may eventually provide the 
evidence for time dilation that has so far eluded detection. 
Subject headings: gamma rays: bursts — galaxies: star formation 



1. INTRODUCTION 

Gamma-ray bursts are among the most distance ob- 
jects observed in the Universe. Having b een detected as 
nearb y as z = 0.0085 for GRB 980425 f Galam a et "ah 
1999t and as dista nt as z = 9.4 for GRB 090429Bl [Cuc° 
chiara et al. 2011) , the redshift range over which these 
objects occur spans almost the entire length of the ob- 
servable Universe. This large dynamic redshift range, 
combined with the transient nature of their prompt 
gamma-ray emission, makes them excellent source with 
which to look for evidence of time dilation due to the 
expansion of the Universe. At such large distances, the 
duration of prompt gamma-ray burst temporal profiles 
should be stretched by a factor of 1 -I- z due to cosmolog- 
ical expansion in Friedman-Lemaitre-Robertson- Walker 
(FLRW) models of the Universe, resulting in an increase 
in the average burst duration as a function of redshift. 

Such evidence for a systematic broadening of GRB du- 
rations as a function of distance has not materialized in 



either the Swift ( [Gehrels et al.| |2004 [) or Fermi dete cted 
GRBs with known redshitt (Campana et al. 2007). In 
fact, the most distance GRBs to date, 
z = 9.4, GRB 090423 at z 



GRB 9804'25 at 
8.1 and GRB 080913A at 



^ 6.7, ha ve measur ed durations of 5.5 sec (Cucchiara 
et al.||2011|), 10.3 sec ( |Salvaterra et al.||2009[ ) and 8.1 sec 
( [Greiner et al.|2009 ) respectively, giving them rest frame 
durations of ~ 0.52 sec, ^ 1.13 sec, and ^ 1.04 sec. 
These values are far shorter than the average rest frame 
duration of ~ 23 — 30 sec for long GRBs detected at 
2 ~ 1. The lack of any correlation between a GRB's du- 
ration and its distance can be seen quite clearly for Swift 
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Fig. 1.— T90 duration vs. redshift for 195 Swift detected GRBs 
as reported via GCN. 

detected GRBs with known redshift shown in Figure [T] 
In this paper, we examine the effects of time dilation 
on the temporal profiles of GRBs through the use of a 
new spectral evolution code that can model the observed 
properties of GRBs as a function of redshift. By us- 
ing prescriptions for the shape and evolution of prompt 
gamma-ray spectra, we can generate a simulated popu- 
lation of single pulsed GRBs at a variety of redshifts and 
observe how their resulting light curves would appear to 
a gamma-ray detector here on Earth. We find that the 
observer frame duration of individual pulses does not in- 
crease as a function of 1 -I- z as one would expect from 
cosmological time dilation. In fact, the duration of indi- 
vidual pulses is seen to decrease as their signal-to-noise 
decreases with increasing redshift, as only the brightest 
portion of a high redshift GRB's light curve is accessible 
to the detector. This effect is analogous to the prob- 
lem faced when measuring the physical radii of galax- 
ies at large distances, where the regions of lower surface 
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Fig. 2. — The hard to soft spectral evolution of i?pij as a function 
of time as prescribed by relativistic kinematics for a spherical shell 
traveling towards the observer with high Lorentz factor. 

brightness fall below the sensitivity of the collecting in- 
strument. In addition, the energy at which most of the 
photons are emitted (i.e. the peak of the z^F^ spectra) 
for high redshift bursts is redshifted to energies where 
traditional gamma-ray detectors become less sensitive, 
decreasing the observed signal even further. The net re- 
sult is that all estimates of duration and energetics for 
high redshift GRBs can only be taken as lower limits to 
their true values. We conclude that GRB temporal pro- 
files that exhibit narrow pulses separated by significant 
quiescent periods, may in fact the tell tail sign of high 
redshift GRBs and predict that the average peak to peak 
duration for a large number of multi-pulsed GRBs as a 
function of redshift may eventually provide the evidence 
for time dilation that has so far eluded detection. 

We present an overview of our population synthesis 
code in Section 2.1, followed by a more in-depth descrip- 
tion in Sections 2.2 through 2.4. We present the result of 
our simulation in Section 3 and discuss the implications 
of our results in Section 4. 

2. GRB MODEL 

Gamma-ray burst continuum spectra can evolve quite 
dramatically over the course of a burst. This evolution 
is generally characterized by an overall softening of the 
spectra, with the peak of the i^F^, spectrum (i?pk) evolv- 
ing through the detector bandpass over the duration of 
the burst. This evolution will be delayed by the effects 
of time dilation for GRBs at high redshifts, resulting in 
a longer observed spectral lag between the high and low 
energy channels and a broadening of the pulse profile. At 
the same time, the observed GRB flux falls as a function 
of increasing luminosity distance. The net effect is that 
soft and faint emission will become increasingly difficult 
to observe with traditional gamma-ray detectors, which 
suffer large drops in sensitivity at low energies. The ob- 
served bursts properties are therefore a complex convo- 
lution of the effects of cosmological redshift and detector 
sensitivity and hence we turn to simulations to obtain 
a better idea of how these bursts would appear in the 
observer frame. 

Two empirical correlations form the basis for the GRB 
model that we have constructed to investigate this ques- 
tion. The first is the hardness-intensity correlation or 
HIC, which relates the instantaneous hardness of the 
spectra and the instantaneous energy flux Fe, within 
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Fig. 3. — The resulting fast rise exponential decay (FRED) pulse 
shape generated by the model by integrating the bolometric GRB 
spectrum at each time bin. 



individual pulses ( Golenetskii et aLllO SSV For the decay 
phase of a pulse, the most common behavior of the HIC is 
a power-law relationship between Fe and the peak of the 
vF^ spectrum, iJpk, of the form Fe oc E^y. , where rj is 
the HIC power-law index. The secon d correlation is the 
hardn ess- fluence correlation or HFC [Liang fc Kargatis| 
( 1996 1 which describes the observation that the instan- 



taneous hardness, or i?pk, of the spectra decays exponen- 
tially as a function of the time- integrated flux, or fluence, 
of the burst. The HFC can be stated as Ep]^ = Eoe-"^/"^", 
where <I>(t) is the photon fluence integrated from the start 
o f the burst and $ p is the exponential decay constant. 

Kocevski et al. ( |2003| have shown that both correla- 
tions can be produced through simple relativistic kine- 
matics when applied to a spherical shell expanding at 
relativistic velocity (i.e. the curvature of a relativistic 
shock front). B y assuming a spectral shape, in this case 
a Band model ^Band et al.|[l9 93), for the GRB's instan- 
taneous photon spectrum and evolving that spectrum in 
time according to these two empirical correlations, we 
are able to reproduce the Fast Rise Exponential Decay 
(FRED) pulse shape that is so ubiquitous in GRB data. 

An example of a time resolved photon spectra and the 
resulting photon flux light curve for a simulated burst 
using a Band model with an initial i?pk,src — 500 keV 
and a fixed a = 1 and f3 = —2.2 can be seen in Figures 
[2] and [3] respectively. The photon flux light curve in Fig- 
ure |3] is calculated by integrating the photon spectra over 
the instrument's energy range and as such is not bolo- 
metric. Therefore the shape of the resulting light curve 
will be greatly effected by the location of E'pk within the 
instrument's energy range. 

In order to convert our modeled photon spectrum into 
a count spectrum and eventually a count light curve, we 
take our simulated photon spectrum and fold it through 
an instrument data response matrix (DRM) . For the pur- 
poses of this analysis, we use a response file fro m the 
Burst and Tran sient Source Experience (BATSE) ( |Mee 



gan et al.|1992| ) that was generated for a real burst which 



occurred nearly at zenith for one of BATSE's Large Area 
Detectors (LAD). The DRM describes the distribution of 
counts over the instrument's energy channels due to the 
arrival of a photon of a given energy. The result is a 
time resolved count spectra as a function of channel en- 
ergy and time. 
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Fig. 4. — The effects of increasing redshift on the observed Ught curve of a GRB pulse. The effects of time dilation are hidden by the 
decreasing signal-to-noise and the pulse duration in fact decreases with increasing redshift. 



We then add a Poisson distributed, energy dependent, 
background spectrum derived from the median back- 
grounds of a sample of BATSE detected bursts to each 
time resolved count spectra. A count light curve with 
a realistic background spectrum can then be produced 
as a function of time by integrating the individual time 
resolved count spectra over the instrument's effective en- 
ergy range for each time bin. The burst duration is 
then measured using a Bayesian block algorithm to de- 
termine periods of emission that are above the Poisson 
background. We refer to this estimate as the TlOO dura- 
tion (as opposed to the traditional T90 duration found 
through flux integration methods). 

We refer the reader to Kocevski et al. (2011) for a 
more comprehensive description of the GRB model used 
for the present analysis. 

3. RESULTS 

Figure H] shows the resulting count light curves for the 
simulatea GRB shown in Figure [2] placed at three dif- 
ferent redshifts, along with the Bayesian block recon- 
struction of the burst duration (solid green lines). Con- 
trary to what is expected due to cosmological time di- 
lation, the duration of the simulated pulse decreases as 
a function of increasing redshift. This is primarily due 
to the decreasing signal-to-noise of the observed pulse 
with increasing redshift, which progressively limits only 
the brightest portion of the GRB's light curve from be- 
ing accessible to the detector. This effect is analogous 
to the problem faced when measuring the physical radii 
of galaxies at large distances, where the regions of lower 
surface brightness fall below the sensitivity of the collect- 
ing instrument. In addition, the energy at which most of 
the photons are emitted (i.e. the peak of the vF^ spectra) 
is redshifted to energies at which the instrument becomes 
less sensitive, reducing the observed signal even further. 

Both of these effects can be seen in Figure [S] where we 
plot the observer frame duration for a set of simulated 
GRB pulses, of equal intrinsic duration but varying in- 
trinsic luminosity, as a function of redshift. The solid 
line represents the expected observer frame duration in 
a FLRW Universe, found by multiplying the pulse's in- 
trinsic duration by (I + 2). The dotted line represents the 
observed duration as measured from the photon flux light 
curves, found by integrating the photon spectra over the 
detector's energy range. The dashed lines represent the 
observed duration as measured from count light curves 
for each of the simulated GRB pulses. These count du- 
rations actually represent the median value of 200 count 



light curve realizations for each redshift bin, computed in 
order to account for statistical fluctuations in the simu- 
lated background which may dominate the duration mea- 
surements as the observed signal-to-noise decreases with 
increasing redshift. The error bars represent the stan- 
dard deviation of the resulting duration distribution for 
each set of 200 realizations. 

The divergence between the predicted observer frame 
duration and the duration of the photon flux light curve 
can be understood as due to the redshifting of the un- 
derlying GRB spectra. Because the simulated photon 
flux light curves are non-bolometric, their shape will be 
greatly effected by the redshifting of towards the 
lower edge of the detector's energy range. Essentially, 
the pulse profile that would have been seen by the detec- 
tor if the burst had occurred at z — 0.01, for example, 
is not the same as the profile seen at z = 2, because the 
lower end of the GRB spectrum is no longer detected by 
the instrument. Therefore, even a perfect detector that 
observes over limited energy range would not faithfully 
measure the expected time dilation effects on a GRB 
pulse as a function of redshift. 

The divergence between the predicted observer frame 
duration and the measured duration becomes wider when 
considering duration measurements made through the 
use of the detector's count data. When noise is added 
to the observed signal, the measured duration turns over 
and begins to decrease with increasing redshift. As the 
pulse's signal-to-noise falls, only the brightest portion 
of pulse becomes accessible to the detector, until the 
observed duration approaches zero and the pulse is no 
longer detected. As seen in Figure [5] the redshift at 
which the transition between a rising and falling dura- 
tion depends largely on the luminosity of the pulse, but 
is also influenced by the burst's intrinsic fi'pk.src, since 
the redshifting of .Bpk.src towards the lower edge of the 
detector's energy range acts to further reduce the signal- 
to-noise of the pulse. 

Transforming the observed durations into the source 
frame by dividing by (1 -I- z) has the effect of further 
widening the difference between the pulse's true intrinsic 
duration and our estimates. The systematic error that 
this duration bias introduces also propagates into our 
estimate of the burst's total energy release, as E^so is es- 
timated by integrating the burst's flux over the observed 
duration. Figures [6] and [7] show both of these effects. 

In FigurelBl we have plotted the estimated source frame 
duration, normalized to the true intrinsic duration, for 
burst's of varying luminosity, but of equal intrinsic du- 
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Fig. 5. — The observer frame duration for a set of simulated 
GRB pulses, of equal intrinsic duration but varying intrinsic lumi- 
nosity, as a function of redshift. Even a perfect detector collecting 
data over a limited energy range would not faithfully measure the 
expected time dilation effects on a GRB pulse as a function of red- 
shift. When noise is added to the observed signal, the measured 
duration turns over and begins to decrease with increasing redshift. 

ration (dashed lines). The turn over observed in Figure 
[5] can now be seen as a steeping of the estimated source 
frame duration as a function of redshift. In ah cases, the 
under-estimation of the source frame duration can ap- 
proach as much as 80% near the detection threshold of 
the instrument. 

The consequences that this systematic error has on the 
estimates of i^iso can be seen in Figure [7] where the ratio 
of the the pulse's estimate to true -Eiso can reach as much 
as 90%, or more, as the pulse nears the detector's sensi- 
tivity threshold. The Eiso measurements presented here 
were k-corrected to a standard energy range of 10 — 10000 
KeV, hence this error in our estimate of E^iso is largely 
due to the under-estimate of the duration over which the 
burst's observed flux is integrated. 

In Figure [Sj we quantify the severity of the duration 
bias as a function of the observed signal-to-noise ratio by 
simulating a 4000 single peak GRBs at z = 1 of equal in- 
trinsic duration, but varying luminosity. The estimated 
observer frame duration normalized to the intrinsic du- 
ration is plotted vs the pulse's observed signal-to-noise, 
with a color coding displaying the peak photon flux of 
the observed signal. Again, the bias in the duration es- 
timate is largest for intrinsically weak bursts with low 
signal-to-noise as seen in the detector. Although the bias 
is far less severe for the highest signal-to-noise bursts, the 
offset between the estimated duration and the true dura- 
tion is still evident and due to the redshifting of -Epk.src • 
The divergence between the true observer frame duration 
and the estimated duration at low signal-to-noise ratios 
reflect the similar breaks seen in Figures [6] and |6] More 
importantly, this divergence occurs at roughly tne same 
signal-to-noise ratio, in this case SNR ^ 25, for all of the 
bursts in this simulation. 



4. DISCUSSION 

The systematic error introduced into our estimate of 
the observer frame duration has severe consequences for 
our subsequent estimates of the source frame duration 
and the isotropic equivalent energy Eiso of a GRB pulse. 
By under-estimating the observer frame duration, any 
transformation into the source frame by dividing by 1 -I- z 
only acts to further diverge our estimates of the source 
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Fig. 6. — The estimated source frame duration, normalized to 
the true intrinsic duration, for a set of simulated burst's of vary- 
ing luminosity. Transforming the under-estimated observer frame 
durations into the source frame has the effect of further widening 
the difference between the pulse's true intrinsic duration and our 
duration estimates 



frame duration from their true values. For the simulated 
pulses shown in Figure |8j this fundamental duration bias 
is largest for bursts detected near the instrument's de- 
tection threshold and can be as much as 90% before the 
pulse is no longer detected by the instrument. Such a 
large divergence between our estimated source frame du- 
ration and the true intrinsic duration result in large sys- 
tematic error in our estimates of the burst's total ener- 
getics. 

Our simulations show that, baring the existence of an 
intrinsic evolution in the duration of GRBs and redshift, 
we do not expect high redshift, single pulsed, GRBs to 
be characterized by extremely long pulse durations. This 
conclusion agrees with the fact that evidence does not ex- 
ist for the systematic broadening of GRB durations as a 
function of redshift in either the Swift or Fermi detected 
GRBs with known redshift. The most distance GRBs to 
date, GRB 980425 at z = 9.4, GRB 090423 at z - 8.1 
and GRB 080913A at z ~ 6.7, have measured durations 



of 5.5 sec ( Cucchiara et al. 2011), 10.3 se c (Salvaterra 
et al.|200 9) and 8.1 sec ( Gremer et al.|2009 ) respectively, 
giving them rest frame durations of ~ 0.52 sec, ~ 1.13 
sec, and ~ 1.04 sec. Our simulations show that these 
values are actually lower limits to the true duration of 
these events and as such, so to are their estimated Eiso 
values. 

Although the results presented in Figure [5] show an 
increase in the observed duration of the five simulated 
pulses at low redshifts, a broad luminosity fu nction, such 
as the one reported by Butler et al. (20101, will act to 
mask such time dilation signatures through the existence 
of intrinsically weak bursts at all redshifts. Our simula- 
tions show that such a broad luminosity function effec- 
tively rules out the possibility that early studies success- 
fully detected signatures of time dilation when consider- 
ing the average duratioi i of BATSE detect ed GRBs as 
a function of peak flux (Norris et al. 1994 ). Th i s is es - 
sentially the same argurnent presented by Band (1994), 
years before the first measurement of redshifts associated 
with GRB afterglows and their host galaxies. 

The fundamental duration bias presented here has 
broad consequences of our understand of the energet- 
ics of GRBs, as the true Ejso distribution may extend 
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Fig. 7. — The burst's estimate normalized to the burst's 

true Eiso for a set of simulated burst's of varying luminosity. The 
systematic error that the duration bias introduces propagates into 
our estimate of the burst's energetics 

to higher energies than had previously bee n estimated 
( [Butler et al.|2010[ [Kocevski fc Butler(2008[ ). It also im- 
plies that attempts to predict the time ot jet tjet breaks 
in afterglow light curve, usually performed by measur- 
ing ii'iso and assuming a canonical coUimation corrected 
energy E^, may also be suspect. If i?iso values for low 
signal-to-noise ratio bursts are truly lower limits to the 
true energy, then their associated tjet values are likewise 
lower limits and their jet break may actually occur at a 
point when the afterglow has faded beyond detectability. 

Although measurements of Eiso may suffer from the 
systematic bias in duration measurements, estimate of 
a burst's peak isotropic equivalent luminosity Liso,pk 
should not be effected, as long as a measurement of 
£'pk,obs allows for a proper k-correction of the burst's 
underlying spectra. This would mean that although Eiso 
evolution studies may in fact be biased at high redshift 
where a greater percentage of detected bursts have low 
signal-to-noise ratios, no such bias should exist when con- 
sidering Liso- Ther efore, studies findin g the lack of lu- 
minosity evolution (Butler et al. 2010) as a function of 
redshift would not be ettected, although the similar at- 
tempts to quantify the evolution of E^so may be suspect. 

Ultimately, it is the shape of GRB pulses that leads 
to the fundamental duration bias observed in this anal- 
ysis. In general, signals with high kurtosis will suffer a 
greater degree of bias in their duration estimate, whereas 
a signals with low kurtosis will suffer no bias at all. For 
example, a duration bias would not manifest for a signal 
exhibiting a step function profile. The duration of the 
step function would simply increase as a function of 1 -I- z 
until it fell below the instrument's detection threshold. 

Although our simulations were performed using a 
BATSE data response matrix, these results should gen- 
erally apply to any flux limited instrument that collects 
time series data. The results presented in Figure [8] indi- 
cate that a burst's observed signal-to-noise ratio can pro- 
vide a useful proxy to determine the severity of the bias 
in the estimates of burst duration. In the simulations 
presented in Figure |8] any bursts detected by BATSE 
with a signal to noise less than ^ 25 should be consid- 
ered heavily biased. In a similar fashion, simulations of 
a large number of GRBs of varying luminosity, duration, 
and source frame spectra for different spacecraft may aid 
in determining the instrument specific signal-to-noise ra- 
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Fig. 8. — The resulting fast rise exponential decay (FRED) pulse 
shape generated by the model by integrating the bolometric GRB 
spectrum at each time bin. 



tio below which duration estimates may be considered 
suspect. 

The existence of the classic bimodal duration distri- 



bution reported in Kouveliotou et al. ( 1993 ) would seem 
to limit the degree to which this duration bias can pro- 
duce arbitrarily short durations for intrinsically long du- 
ration GRBs, despite that we have shown that GRB du- 
rations can be underestimated in some cases as much as 
80—90%. This apparent contradiction can be reconciled 
by considering that most GRBs are in fact not single 
pulsed events, but comprise multiple emission periods. 
Contrary to the effects that increasing redshift have on 
the observed duration of individual pulses, the duration 
between pulses, i.e. the quiescent periods, should in fact 
increase. Therefore, although long GRBs that exhibit a 
single, low signal-to-noise, pulse may result in observed 
durations that would allow them to be confused with the 
short GRB population, a majority of long GRBs would 
not suffer this fate even though their durations would 
still be under-estimated. None-the-less, our work demon- 
strates that attempts at distinguishing between long and 
short GRBs based on a burst's tempo ral properties alone 
( [Fynbo et al.||2006[ [Ofek et al.|[2007l )) cannot be trusted 
and additional information about tne burst (e.g. spec- 
tral lag, host galaxy properties, hardness ratios) must be 
brought to bare. 

A profile exhibiting narrow pulses separated by signif- 
icant quiescent periods, may in fact be the tell tail sign 
of a high redshift GRB. Although there are undoubtably 
low redshift GRBs which intrinsically exhibit such tem- 
poral profiles, additional spectral information could be 
used in order to distinguish these bursts from their high 
redshift counterparts. Therefore, a burst's hardness as 
measured by an instrument with broad-band energy cov- 
erage, such as the Fermi-GBM, combined with its ob- 
served temporal profile and signal-to-noise ratio could 
play an important role in distinguishing high redshift 
events. 

Finally, we conclude by predicting that the average 
peak-to-peak time for a large number of multi-pulsed 
GRBs as a function of redshift may eventually provide 
the evidence for time dilation that has so far eluded de- 
tection. 

5. CONCLUSIONS 
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By simulating the effects of time dilation on the evo- 
lution of the prompt gamma-ray spectra, wc have found 
that: 

• The observer frame duration of individual GRB 
pulses does not always increase as a function of redshift 
as one would expect from the cosmological expansion of 
a Priedman-Lemaitre-Robertson- Walker Universe. 

• The observed duration of individual pulses can 
actually decrease with increasing redshift as only the 
brightest portion of a high redshift GRB's light curve 
become accessible to the detector. 

• The net result of this fundamental duration bias is 
that all estimates of duration and energetics for high 
redshift GRBs can only be taken as lower limits to their 
true values. 

• Bursts detected near the instrument's detection 
threshold are effected the most by this bias and an 
instrument specific signal-to-noise ratio below which 
duration estimates may be considered suspect may be 
isolated through simulations. 

• Any attempts at distinguishing between long and 
short GRBs based on a burst's temporal properties 



alone should not be trusted. Additional information 
about the burst its spectral lag, host galaxy properties, 
and hardness ratios must be used in such classification 
attempts. 

We conclude that GRB temporal profiles that exhibit 
narrow pulses separated by significant quiescent periods, 
may in fact the tell tail sign of high redshift GRBs and 
predict that the average pcak-to-pcak duration for a large 
number of multi-pulsed GRBs as a function of redshift 
may eventually provide the evidence for time dilation 
that has so far eluded detection. 
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